Landau theory of the Fermi-liquid to electron glass transition 
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A lattice model of spinless interacting electrons is used to formulate the Landau theory of the 
Fermi liquid to electron glass quantum phase transition. We demonstrate that the presence of 
additional random site energies does not affect the character of the transition, once the replica 
symmetry breaking is considered self-consistently at the mean-field level. Inside the glass phase, 
the low temperature conductivity assumes a non-Fermi liquid 8a ~ y 3 / 2 form, in agreement with 
recent experiments. 

PACS numbers: 75.20.Hr, 71.55.Jv 
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Numerous recent experiments have demonstrated that 
the phase diagram of the low carrier density systems 
proves remarkably rich even in the absence of magnetic 
fields Q| . The fascinating strong correlation physics has 
been observed whenever the energy scale of the Coulomb 
repulsive interactions is comparable to that of the Fermi 
energy. Once disorder is weak enough, the interactions 
alone can lead to insulating behavior, characterized by 
hard gaps in the charge sector. But if the magnitude 
of random potential is comparable to the Fermi energy, 
electrons approach the strongly localized regime. In this 
case, the interplay between localization and Coulomb in- 
teractions is generally expected to lead to (gapless) glassy 
ordering of electrons |§, [| . 

Very interesting new information on this problem 
comes from recent experiments by Bogdanovich and 
Popovic Q. In this work, strong long-time fluctuations 
of conductivity were observed as the electron density is 
reduced below some critical value, indicating a dramatic 
slowing down of electron dynamics. This behavior was 
attributed to the presence of the glassy freezing that ap- 
pears to come from the charge degrees of freedom. Inter- 
estingly, the glassy behavior seems to emerge appreciably 
before the metal-insulator transition, thus identifying an 
intermediate metallic glass phase. Inside the glass phase, 
unusual temperature dependence of the conductivity was 
observed, which was fitted to the form 5a ~ T 3 / 2 . 

The main goal of the current Communication is to the- 
oretically examine the nature of the quantum phase tran- 
sition from a Fermi liquid to such a metallic glass phase, 
and the account for the resulting modifications of charge 
transport. Previous works have extensively examined the 
quantum paramagnet to spin-glass transitions, regarding 
thus the spin degrees of freedom as the only relevant ones 
at criticality|| Q, @, §|. Here, we use a similar approach 
to study the onset of glassiness in the charge sector. Us- 
ing a recently developed dynamical mean-field formula- 
tion ||, we construct a Landau theory which provides a 
complete description of the electron dynamics near the 
relevant quantum critical point. 

We consider a lattice model of spinless interacting elec- 
trons at half-filling in the presence of on-site randomness, 



as given by the Hamiltonian 



ijCjCiCjCj. 



(1) 



Here t%j denotes the corresponding hopping elements, 
while Vij describes the inter-site Coulomb interactions. 
The distribution of the random site energies e% is as- 
sumed to be Gaussian, with variance W. As the simplest 
model H for the glassy freezing of electrons, we choose 
the inter-site interactions also to be Gaussian dis- 
tributed random [^| variables, with variance V . For this 
model, it is straightforward || to employ the methods 
of dynamical mean-field theory |l0| which is formally ex- 
act in the limit of large coordination. The averages over 
randomness are carried out with using standard replica 
methods |ll| , which are also used to identify the emer- 
gence of the glassy phase. This procedure leads || to the 
following single site effective action (w m = 2ir(m + 1)T) 

Ses(i) = y^y^ ct a (q; m ) [iuj m + t 2 G(uj m )] cf(cj m ) 

- f f dr x dr 2 J ^ V 6n%(n)Q ab (n " r 2 )5n\{T 2 ) 
Jo Jo { 2 ab 

+ y^K(r 1 «(r 2 ) (2) 

ab ) 

In the formula above, the functional integration is per- 
formed over the fermionic fields c"(r) (a — 1, ...,n are 
the standard replica indices). Snf(r) — cl a (T)cf(r) — | 
denotes the deviation of the density from half-filling. We 
took into account that the "Weiss" (cavity) field has the 
form W s (t\ — t 2 ) = t 2 G(r\ — t 2 ) for electrons on a Bethe 
lattice. The Green function G{t\ — t 2 ) along with the 
order parameter Q ab {ri — t 2 ) must be determined self- 
consistently using the effective action given by Eq. (||) || . 

G(r 1 -r 2 )=<4 a (TiK(r 2 )>eff, (3) 
Q ab {ri - 75) =< <5<(n)<5nf (r 2 ) > off • (4) 

The replica diagonal components of Q aa (T) represent the 
averaged dynamic compressibility, while the parameters 



2 



Q (a ^ b) are time independent and related to the fa- 
miliar Edwards- Anderson (EA) order parameter^). Sim- 
ple analysis of Eqs. shows that q a b is non-zero 
everywhere, once W ^ 0. This is the consequence of the 
non-uniform density due to the on-site random potential. 
As a result, the glass transition for W ^ assumes the 
character of a De Almeida- Thouless line M, |ll), where 
special care is needed in formulating a Landau theory. 
As it will be clear from below, the relevant Landau the- 
ory can be rigorously formulated, only if W is small, but 
our conclusions remain qualitatively correct for arbitrar- 
ily large W. If V is much smaller than t, action Eq. (|^) 
describes the usual disordered Fermi liquid fl3| , while in 
the opposite limit the glassy ordering persists down to 
T = 0||. The quantum phase transition to the glassy 
phase occurs in this model at some critical value (t/V) CT , 
which has weak dependence on W ||. 

To obtain the Landau functional we must perform a 
cumulant expansion in Eq. (|J) treating the term with 
Q ab (ji — r 2 ) as perturbation. Before doing this, it is 
necessary to shift the Q-matrices, eliminating thus the 
non-critical regular part in diagonal elements ((3 = 1 /T) : 



replica symmetry breaking. To make this point clear we 
first identify the part of the action that contains only q a f,: 

ft = -R^rq 2 ^Trq 3 - ?± ^( 9ab ) 4 - R A £ g o6> (8) 



where, 



Ri = (3(D(Q) + /3<7ea), i? 2 = /3 2 



i?4 



(3w 2 v 2 



(9) 



Fermi liquid phase. As we emphasized previously, 
there is no replica symmetry breaking in the Fermi liquid 
phase. Therefore, it is natural to choose the parametriza- 
tion q a b = <7ea in this phase and take the variational 
derivative of Eq. (M) with respect to <7ea, obtaining thus: 



2v W 2 V 2 
2D(0)q EA + ^Lq 3 EA + — ^— = 



(10) 



Varying subsequently Eq. (g) with respect to D(uj n ) we 
arrive at the equation: 



Q ab (co n ) -> Q afc K) - JC5 ab 



W 2 



(5) 



r + \uj n \ + u 



u n ) + qEA 



D 2 {uj n ) 



Constant /C must be formally determined from the con- 
dition of absence of the term XL„ SnJ ( 3 ab (<- t; «)] 2 i n the 
underlying Ginzburg-Landau actionpf, that reads: 



/^=E 



r + \u) n \ 
V 2 



Q aa (tO r , 
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EE^ ai, ( w n)Q 6c K)Q ca (o 



E^ a >«) 

Q a6 K = 



v 4 

2y 1 

' W 4 (3 2 



E D(u 1 )D(lj2)D(-l} 1 - wa - w n ) =(0,1) 



abc uj n 



that closes the system of equations determining qEA and 
J?(uj n ). All dangerous terms proportional to (3 vanish 
'because we have judiciously chosen q a b equal to <7ea from 



dTidT 2 J2[Q ab (n-T 2 )]\ 



a b 



Here r, being some function of t/V, is the parameter 
that governs the transition, while u and y are taken at 
(t/V) cr . The presence of the last term, responsible for 
the RSB instability, is crucial to further analysis. Ac- 
cordingly, we employ the following mean-field ansatz for 
the Q-matrices: 



the very beginning. Note that Eq. (|llj) does not contain 
(6)W at all, while qea = only when W = 0, as can be 
seen from Eq. ( |l0| ) 

Electron glass phase. The saddle point solution of Eq. 
(||) in this phase must be characterized by the Parisi func- 
tion q(s) with < s < 1. Considerations, that are com- 
pletely analogous to the classical case §, m n> iead to 

the functional form that has two plateaus: 



V' z Q ab (cu n ) 



D(tu n ) + /3<?ea<L„,o, a = 6, 
/3<7a&<L„,o, a^b. 



(7) 



In Eq. (^) <7ea is the EA order parameter, and it is 
assumed that q aa = 0. The (3- prefactors are chosen 
to ensure the finite limit of the free energy density as 
T — > 0. We must insert Eq. (0) into the action Eq. (||) 
and find the saddle point solution with respect to the 
variations of q a b, qEA and D(u) n ). q^A and q a b should 
not, however, be varied independently 12 . They must 
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(12) 



R 2 - (R 2 - 4 J R 1 fl 3 ) 1 / 2 ^ . 
Qea = 7T5 > s i < s < 1. 



obey an additional relation between them that depends 
on the presence of the glassy ordering and, hence, the 



The function q(s) in the above equation saturates at the 
value q(l) = qEA- This parametrization is in agreement 
with the definition 5ea = m aXa/&9a&- Combining Eqs. 
(0) and (11^) we obtain that in the glassy phase 



Qea 



-[D(W 4 }/y, 



(13) 
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This equation connecting D(0) and <?ea plays the same 
role as Eq. ( |lC)| ) in the Fermi liquid phase. We must 
substitute then the solution Eq. ( p^ ) into Eq. (@) and 
add subsequently the result to the remaining part of the 
action that contains D(u> n ). We will not write down 
the expression for the free energy density obtained after 
lengthy, but straightforward calculations, that are largely 
identical to those done in the Appendix C of Ref. p2[ . 
We will state only that as a result of condition Eq. ( |l3| ) , 
there will be no terms in the free energy in which D(Q) 
and qEA are coupled explicitly. Taking the variational 
derivative with respect to D(ui n ) results in the same Eq. 
(|i"T|), albeit in the glassy phase qea is connected with 
D(0) by means of Eq. (|l|) rather than Eq. ©. 

Approaching the glass transition. The necessity to ob- 
tain -D(O) from Eq. ( p"l| ) leaves us with the task of its 
self-consistent solution. The exact analytical solution of 
this non-linear integral equation is clearly out of ques- 
tion. However, close to the T — transition point the 
leading order of the correct solution is possible to obtain. 
The approximations we use hinge also upon the smallness 
of W and, consequently, q^A- 

We notice first that, if y = 0, the complete solution is 
well-known to be§, |, § D(ui n ) = -^jc^fTA, with A 
turning to zero right at the critical point. Let's assume 
that for y ^ the leading approximation of D(uj n ) con- 
tains the same square root singularity as for y = 0, and 
analyze the role of the last two terms in Eq. (|ll]) . The 
key point in this analysis is the value of the integral 



J(fi m ) =TJ2 y/\u n \ + Ay/\u n + n m \+A, (14) 

that at T = is simply calculated to be: 

-.2 oa r>2 IOJ + 2A 



i ,,,, n 2 2A W n 2 

J (ST) = — I A 3 / 2 In —4 + — Arch- 



101 



101 



A arcsin ■ 



mi 



01 +2A 



(A) 3/2 ^ + A (15) 



In Eq. (|T^) A w denotes the upper critical cutoff of the 
order of unity. We see that all the terms in J(O), ex- 
cept the first one proportional to A 3 / 2 , are of the or- 
der of 0(A 2 ,S1 2 ). This means that the prelast term in 
Eq. ([ll]) gives contributions that depend quadratically 
on small parameters A and u> n , and, thus, subdominant 
to the leading terms that scale linearly with them. The 
cutoff-dependent part in its turn leads to a mere renor- 
malization of the coefficient u in uqEA- We denote this 
renormalized term as uqEA- Inserting then J{lj\ + uj n ) 
to the last term of Eq. (^lj) and integrating over L02, we 
conclude, that the only effect this term produces is to 
renormalize the critical value r c . This allows us to omit 
formally the last two integrals in Eq. (|ll]) and resolve the 
ensuing quadratic equation for D(o; Tl ), obtaining that in 
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FIG. 1: Zero-temperature phase diagram demarcating the re- 
gions of Fermi liquid and electron glass phases. For small W , 
r — r c the curve separating these two phases scales roughly 
as W ~ (r c — r) 3,/2 . The limiting behavior of the EA order 
parameter qba and gap A in different regimes is presented in 
the text. 



the Fermi liquid phase: 



D(w n ) = - v / kjTA, 



uqEA 



(16) 



(17) 



in the leading approximation. It is easily verifiable d pos- 
teriori that the first term in Eq. ( )l6| ) , being inserted into 
the last two integrals of Eq. (]Tl]) , renders the contribu- 
tions of the higher order of smallness compared to uqea- 
Together Eqs. (|o|) and ( |l6|) determine completely the 
T = behavior of the disordered phase near the quan- 
tum critical point. As as a result of their solution, one 
can distinguish the following regimes on a (r — r c , W) 
plane, schematically depicted on Fig. (1). 

(I) In this regime, in which W <§; (r — r c ) 3 / 4 and can 
be treated as a perturbation, (/ea = {W 2 V 2 )/A^fr — r c , 
A « r — r c . 

(II) This region is characterized by |r — r c | 3 / 4 <§; W . As 
a result, we have <7ea « {W 2 V 2 /4^/u) 2 ^ 3 and A « uqEA 
in this regime. 

(III) This regime, in which (r c — r) 3 / 4 ^ W, is the 
closest to the T = critical boundary. EA order param- 
eter, that crosses over to its value in the glassy phase, is 
given by q EA = [(r c - r)/u] + (A/u), with 



2y(r c 



M 2 V 4 



4(r c 



(18) 



From Eq. ([18]) it is easily seen that A vanishes at the 

critical line given by W = (8y/3) [(r e - r)/uV 2 f' 2 . This 
is exactly the line that determines the transition to the 
glassy phase obtained from the joint solution of Eqs. 
@, @ and (fD|. Taking into account Eq. (|l3|), we 
resolve similarly the equation for D(u) n ) to get that be- 
low this line 



D(u n ) 



v 4 



qEA 



(19) 
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Within the spin glass phase, that obtains for all values of 
W , A is zero everywhere. Therefore, the imaginary part 
of the local dynamic susceptibility has the non- Fermi liq- 
uid singularity of the form ~ sgn(o;)w 1 / 2 . The transition 
to the electron glass phase occurs as a second-order tran- 
sition and is of the same character for both zero and 
non-zero on-site randomness W . 

Finite temperature behavior. The evaluation of the 
temperature-dependent correction to the integral J(f2 m ) 
in Eq. (|l4j) leads to the following results. If A <C T, the 
correction scales oc T 2 for \Q m \ ~ T and oc T 3 / 2 |f2 m | 
for |n ro | > T. If A > T, it behaves cx T 2 for 
|fi m | ~ T and cx T 2 \fl m \/y/A for |fi m | > T. As a 
result, the equations governing the dependence of pa- 
rameters qEA(T) and A(T) at finite temperature can be 
obtained from those considered above by the formal sub- 
stitution of r c by r c (T) = r c - c x T 2 j\J~K for A > T, and 
r c (T) — r c — C2T 3 / 2 in the opposite limit. c\ and C2 are 
now, however, some complicated cutoff-dependent coef- 
ficients of the order of unity. Those modifications lead 
to the myriad of limiting cases close the quantum criti- 
cal point that will not be enumerated here. Instead, we 
mention that only if T is small enough in the Fermi liq- 
uid phase, all temperature corrections behave as T 2 with 
some large prefactors arising because W, \r — r c \ <C 1. 
On the contrary, for all temperatures in the glassy phase 
q EA (T) = [r c -r - c 2 T 3 / 2 ]/u. The calculation of the 
electron self-energy Q] in the lowest order using Eqs. (||), 
(0) and suggests, that the leading temperature cor- 
rection to the elastic scattering rate is proportional to 
T 3 / 2 in the glassy phase. This gives rise to the non- 
Fermi liquid temperature dependence of the conductivity, 
5o~ ~ T 3 / 2 , in qualitative agreement with the experimen- 
tal observations [0. 

We emphasize that our results are based on the non- 
perturbative treatment of the RSB term in Eq. (|S ) and 
the specific conditions, connecting q^A and q a b 12 . An 
alternative saddle point solution of Eq. (||) , in which <7ea 
and q a b are allowed to be varied independently || , leads 
the gapless glassy phase, only if W = 0, which we be- 
lieve is incorrect. Our solution which has zero gap for 
all finite W, seems to be the one in agreement with the 
recent numerical simulations fl6| . Though our theory is 
formulated for small W, the saddle-point Eqs. ([ll]) and 
( |i~3| ) are free from this parameter. This strongly suggests 
that the similar set of equations, with renormalized co- 
efficients, should describe the glassy phase close to the 
transition for all W, including the limit W — * 00 ||. 

To conclude, we have presented a Landau theory de- 
scription of the disordered Fermi liquid to electron glass 
quantum critical behavior. Our results represent an ex- 
act solution frTsfl of the model within a dynamical mean- 
field formulation, which is formally exact in the limit of 
large coordination. It is important to note that a glass 
transition having a character of a De Almeida-Thouless 
line plj, such as the one we describe, generally emerges 



within mean-field models. An alternative formulation 
Q, based on droplet approaches predicts the absence 
of such transitions for models with short-range interac- 
tions. In the case of an electron glass, the existence of 
the long-range Coulomb interaction opens a possibility 
that droplet approaches are not relevant, and that the 
glassy behavior of electrons could be well described using 
mean-field models. This possibility seems to find support 
in very recent experiments p0[ , which provide striking 
evidence of scale-invariant dynamical correlations inside 
the glass phase, consistent with the hierarchical picture 
of glassy dynamics as emerging from mean-field models. 
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